Superconductors with Magnetic Impurities: Instantons and Sub-gap States 
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When subject to a weak magnetic impurity potential, the order parameter and quasi-particle 
energy gap of a bulk singlet superconductor are suppressed. According to the conventional mean- 
field theory of Abrikosov and Gor'kov, the integrity of the energy gap is maintained up to a critical 
concentration of magnetic impurities. In this paper, a field theoretic approach is developed to 
critically analyze the validity of the mean field theory. Using the supersymmetry technique we find 
a spatially homogeneous saddle-point that reproduces the Abrikosov-Gor'kov theory, and identify 
instanton contributions to the density of states that render the quasi-particle energy gap soft at any 
non-zero magnetic impurity concentration. The sub-gap states are associated with supersymmetry 
broken field configurations of the action. An analysis of fluctuations around these configurations 
shows how the underlying supersymmetry of the action is restored by zero modes. An estimate of 
the density of states is given for all dimensionalities. To illustrate the universality of the present 
scheme we apply the same method to study 'gap fluctuations' in a normal quantum dot coupled 
to a superconducting terminal. Using the same instanton approach, we recover the universal result 
recently proposed by Vavilov et al. Finally, we emphasize the universality of the present scheme for 
the description of gap fluctuations in d-dimensional superconducting/normal structures. 



PACS numbers: 74.62.Dh, 71.55.-i, 74.40.+k 



I. INTRODUCTION 



While the spectral properties of a singlet s-wave super- 
conductor are largely unaffected by weak non-magnetic 
impurities 0, the pair- breaking effect of magnetic im- 
purities leads to the gradual destruction of superconduc- 
tivity. Remarkably, the suppression of the quasi-particle 
energy gap is more rapid than that of the superconduct- 
ing order parameter, admitting the existence of a narrow 
'gapless' superconducting phase in which the quasi- 
particle energy gap is destroyed while the superconduct- 
ing order parameter remains non-zero. Now, according 
to the conventional (mean-field) description formulated 
in the seminal work of Abrikosov and Gor'kov (AG), an 
energy gap is maintained up to a critical concentration 
of magnetic impurities (at T = 0, 91% of the critical con- 
centration at which superconductivity is destroyed). Yet, 
being unprotected by the Anderson theorem, it seems 
likely that the gap structure predicted by the mean-field 
theory is untenable and must be destroyed by 'optimal' 
fluctuations of the random impurity potential. Indeed, 
since the pioneering work of AG, several authors |^-||] 
have explored the nature of 'sub-gap' states in the super- 
conducting system. The aim of this work is to present a 
detailed investigation of the spectrum and profile of sub- 
gap states in superconductors subject to a weak mag- 
netic and non-magnetic impurity potential, thus system- 
atically improving upon the mean-field theory of AG. Our 
preliminary findings have already been reported in a re- 
cent letter ||. 

In the earliest works on the subject ]3|-||, attention 
was focussed on the the influence of strong magnetic 
impurities. In particular, in the unitarity limit, it was 
shown that a single magnetic impurity leads to the local 



suppression of the order parameter and creates a bound 
sub-gap quasi-particle state (^] . For a finite impurity con- 
centration, these intra-gap states broaden into a band 0] 
merging smoothly with the continuum bulk states. 

By contrast, starting with a weak magnetic impurity 
distribution (i.e. one in which the magnetic scatter- 
ing can be treated within the Born approximation), the 
mean-field theory of AG Q predicts a gradual suppres- 
sion of the quasi-particle energy gap. Defining the di- 
mensionless parameter 
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where |A| represents the value of the homogeneous self- 
consistent order parameter, and r s denotes the Born scat- 
tering time due to magnetic impurities, the AG theory 
shows the gap to follow the relation 



£ gap (r s ) = |A| (l - C 
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showing an onset of the gapless region at £ = 1 (note 
% = 1 throughout). Within the same mean-field theory, 
for £ < 1, the self-consistent order parameter varies as 
|A| = | A| exp[— 7r£/4] where |A| represents the order pa- 
rameter of the clean superconductor, confirming that the 
order parameter is finite at the onset of the gapless phase. 
The precise variation of the quasi-particle energy gap is 
compared to that of the self-consistent order parameter in 
Fig. |l|. Staying within the framework of the mean-field 
theory, one can obtain a smooth interpolation between 
the strong and weak impurity scattering behaviors ||] . 
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FIG. 1. Variation of the energy gap -E gap and the 
self-consistent order parameter |A| as a function of (normal- 
ized) scattering rate 2/t s |A|. 

However, even for weak disorder, it is apparent that 
optimal fluctuations of the random potential must gen- 
erate sub-gap states in the interval < £ < 1, and 
therefore provide non-perturbative corrections to the self- 
consistent Born approximation used by AG. Extending 
the arguments of Balatsky and Trugman J7|, a fluctua- 
tion of the random potential which leads to an effective 
Born scattering rate 1/t' s in excess of 1/t s over a range 
set by the superconducting coherence length, 
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induces quasi-particle states down to energies E e&v (t' s ). 
These sub-gap states are localized, being bound to the 
region or 'droplet' where the scattering rate is large. 

The situation bears comparison with band tail states 
in semi-conductors. Here rare or optimal configurations 
of the random impurity potential generate bound states, 
known as Lifshitz tail states which extend below 
the band edge. However, the correspondence is, to some 
extent, superficial: band tail states in semi-conductors 
are typically associated with smoothly varying, nodeless 
wavefunctions. By contrast, the tail states below the 
superconducting gap involve the superposition of states 
around the Fermi level. As such, one expects these states 
to be rapidly oscillating on the scale of the Fermi wave- 
length A f , but modulated by an envelope which is local- 
ized on the scale of the coherence length £. This differ- 
ence is not incidental. Firstly, unlike the semi-conductor, 
one expects the energy dependence of the density of 
states in the tail region below the mean-field gap edge 
to be 'universal', independent of the nature of the weak 
impurity distribution, but dependent only on the pair- 
breaking parameter £. Secondly, as we will see, one can 
not expect a straightforward extension of existing theo- 
ries |l(]|n[] of the Lifshitz tails to describe the profile of 
tail states in the superconductor. 



In the BCS approximation, the random system we con- 
sider is specified by the Gor'kov Hamiltonian 
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where the matrix components index the particle/hole 
content, and 
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e F + V(r) + JS(r) ■ <j s 



denotes the normal component of the random Hamilto- 
nian. Here we take A = ^(c^Cf) to be spatially homoge- 
neous and determined self-consistently from the conven- 
tional AG mean-field theory. We will assess the validity 
of this assumption below. In addition to the weak poten- 
tial impurity distribution V(r), the particles experience a 
quenched random magnetic impurity distribution JS(r) 
where J represents the exchange coupling and Pauli ma- 
trices {crf p } operate on the spin indices. The magnetic 
S(r) and non-magnetic V(r) random impurity potentials 
are both taken to be Gaussian (^-correlated with zero 
mean and variance 



(JS ct (r)JS /3 (r')> 
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respectively, and v represents the average density of 
states (DoS) per spin of the normal system. 

To keep our discussion simple, and to make contact 
with the AG theory, we will take the quenched dis- 
tribution of magnetic impurities to be 'classical' and 
non-interacting throughout. For practical purposes, 
this entails the consideration of structures where both 
the Kondo temperature p] and, more significantly, the 
RKKY induced spin glass temperature [O are smaller 
than the relevant energy scales of the superconductor. 
The remaining energy scales are arranged in the quasi- 
classical and dirty limits: 



e F »l/r»(|A|,l/r s ) 
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where r represents the transport time associated with 
non-magnetic impurities. We remark that these limits 
are not compatible with the situation in which the mag- 
netic impurities provide the only source of scattering, 
i.e. £ ~ I = vpT s . We should not, therefore, expect a 
straightforward comparison with the analysis of Ref. M . 

Within the approximations above we will derive a 
quasi-classical field theory of the disordered supercon- 
ducting system. Following the approach of Ref. we 
will express spectral properties of the system in terms of 
an intermediate energy scale action which accommodates 
the quantum interference properties of the superconduct- 
ing system. By investigating stationary inhomogeneous 
instanton field configurations of the action, we will ex- 
pose the structure and profile of the sub-gap states, and 
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thereby obtain an analytical expression for the DoS to 
exponential accuracy. Finally, we will comment on the 
universality of the present scheme by applying the same 
technique to investigate 'gap fluctuations' in a quantum 
dot coupled by open channels to a superconducting ter- 
minal. 

Although the analysis is straightforward, the technol- 
ogy is somewhat involved. We have therefore decided to 
summarize the main conclusions of this investigation here 

I 



i/(e < E g 



exp 



in the introduction. In particular, we will show that, in 
the vicinity of the mean-field gap edge, the sub-gap DoS 
of the (/-dimensional system is dominated by tail states 
which are confined to droplets of size 
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The corresponding sub-gap DoS takes the form 
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where g = vDL d 2 denotes the bare dimensionless con- 
ductance and f d (0 = a d C- 2/3 (l - C 2/3 )" (2+d)/s repre- 
sents a dimensionless function of the control parameter £ 
(ad const.). When reparameterized in terms of the DoS 
just above the mean-field gap edge Q 
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where 



a; 3 / 2 = 4 ™i<y ^ c - 2/3 (i " c2/3yl/i 

the expression for the sub-gap DoS can be brought to the 
more compact form 
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with ad some numerical constant. 

In the zero-dimensional system, although the interpre- 
tation of the optimal fluctuation as a localized droplet 
is no longer appropriate, the expression above correctly 
interpolates to d = and coincides with the universal ex- 
pression for gap fluctuations proposed in Ref. The 
surprising dependence of the result on the dimensionless 
distance from the mean-field gap is one of the reasons 
why a Lifshitz argument appears difficult to construct 
for this problem. 

The paper is organized as follows: in section [n| a theory 
of the statistical properties of the Gor'kov Green function 
is developed within the framework of a supersymmetric 
field theory involving a non-linear cr-modcl functional. 
Here we follow closely the analysis of Ref. |l3| and 
(c.f. Ref. 0)). As a result we identify the conventional 
AG mean-field equation with the homogeneous saddle- 
point equation of the effective action. In section [II we 



show that the non- vanishing of the DoS beneath the gap 
predicted by the standard AG theory is due to the ap- 
pearance of inhomogeneous "instanton" saddle points of 
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finite action. We identify the profile of these instantons 
with the envelope modulating the quasi-classical sub-gap 
states. The instanton configurations considered break 
the underlying supersymmetry of the action, and a de- 
tailed examination of fluctuations is required to under- 
stand how supersymmetry is restored by zero modes, as 
well as to appreciate fully how these configurations are 
able to contribute to the DoS. In section IV we examine 



the zero dimensional limit of the problem and compare 
our results to the recent literature [g 15|. In doing so, we 
provide an explanation of the universal results reported 
in Ref. p5[ , and discuss the universality of the d > 
result. Finally, in section [v] we speculate on potential 
generalizations of the results presented here. 



II. FIELD THEORY OF THE DISORDERED 
SUPERCONDUCTOR 

The construction of the field theory of the disor- 
dered superconductor follows the quasi-classical method 
of Eilenberger JlSj and Usadel |nj elevated to the level of 
an effective action. The starting point of the analysis is 
the generating functional for the single-particle Gor'kov 
Green function for the non-interacting quasi-particle Bo- 
goluibov Hamiltonian. Here we borrow our notation from 
Ref. [ffl. 



A. Generating Functional 

Single-particle properties of the Gor'kov Hamilto- 
nian (0) are obtained from the generating functional 



Z[J] 
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where e_ = e — iO and the supervector fields have 
the internal structure i/j = ( i/>j tpi V| ) > '0 T — 
( ipl V*t "'Pi)- As usual, by introducing both com- 
muting and anticommuting elements the normalization 
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2[0] = 1 is assured. (The generalization to the con- 
sideration of higher point functions, which involves the 
extension of the field space, follows straightforwardly.) 

To condense the notation, it is convenient to perform 
the rotation ip i— ► ip' = Utp, ip i— ► r/r = ipU' with 



U = 



1 
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after which the Gor'kov Hamiltonian takes the form 



H = I + V(r) 
1 2m w 



e F ) ® v7 + JS(r) • ct s 



The unusual phase coherence properties of the super- 
conducting system rely on the particle/hole or charge 
conjugation symmetry 
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To easily accommodate these effects, it is convenient to 
introduce the further space doubling, 



2ip(H - e_)^ = ip(H - e-)ip + ip 1 (H 1 - e_)^ J 
= ${H - e_)V> + ^ t (ct™ ® of Hog ® <J™ ~ t-W 
= 2%>(H- e-(T3 C )* (9) 



where, defining the Pauli matrix a^ c which operates in 
the charge-conjugation (cc) space, 
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This completes the formulation of the generating func- 
tional as a field integral involving 16-component super- 
vector fields $ and The latter obey the symmetry 
relations 
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(10) 



with 7 = E BB icr% c — i? FF crf c , where Ebb — diag(l,0) BF 
and E FF = diag(0, 1) BF project into the Boson-Boson and 
Fermion-Fermion sectors respectively. 



B. Ensemble Averaging 



Cast as a field integral, the impurity average of the gen- 
erating functional over the Gaussian distributed random 
impurity potentials is straightforward. Separating the 
regular from the disordered components of the Hamilto- 
nian H = Ho + V(r)a^ H + JS(r) -<j sp , an ensemble average 
over the random potentials obtains 
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The interactions generating by the impurity averag- 
ing can be decoupled by the introduction of a Hubbard- 
Stratonovich field. Beginning with the non-magnetic dis- 
order, slow modes of the action are identified by rewriting 
the action in the approximate form 
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AlTVT 



dr 
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2'kvt 



]T str C(-q)C(q), 



|q|<£ 



where <(q) = ^ k *( k ) ® *(- k + <ti a T ■ The latter can 
be decoupled by the slowly varying 16 x 16 supermatrix 
field Q(r) according to the identity 



exp 
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= J DQ exp ^ s tr^Q(q)Q(-q)-i-Q(q)C(-q) 



The symmetry properties of the fields Q(r) are inherited 
from the dyadic product \I/(r) ^(i^cr™. Making use of 
Eq. ( |Io| ) one finds the symmetry relations 



Q = a\ n ® erf 7<3 T 7"V PH <g> a. 



(12) 



The interaction generated by the magnetic impurity 
averaging can be treated |^0| by performing all possible 
pairings and making use of the saddle-point approxima- 
tion Q(r) = 2(f(r) ® *(r)of H )<p/7T2/. This leads to the 
replacement 



1 
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Such an approximation, which neglects pairings at non- 
coincident points is allowed by the strong inequality 
(£/£) d <C 1. In addition we discard the contraction 
(\jr<j sp \]/^. xhe term generated by this procedure could 
in any case be decoupled by a slow Bosonic field S(r) 
which would immediately be set to zero for the singlet 
saddle-points that will be the basis of this paper. 

Gaussian in the fields \& and "J, the functional integra- 
tion can be performed explicitly after which one obtains 
(Z[0]) v ,s = J DQexp(-S[Q}) where 
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Further progress is possible only within a saddle-point 
approximation. Following Ref. p3[ , the saddle-point 
analysis will be conducted in a two-step process: 



C. Saddle-point Approximation and the Non-linear 
(7-model 

The first task is to make use of the quasi-classical pa- 
rameter tpT 3> 1 to construct an intermediate energy 
scale action. Dropping the symmetry breaking perturba- 
tions e, |A| and 1/t s , a variation of the action with re- 
spect to fluctuations of Q obtains the saddle-point equa- 
tion 



Q(r) = — (r 



+ h + iQ/2T]- 1 \r) 



where kg = p 2 /2m — e F . Taking into account the ana- 
lytical properties of the Green function and the quasi- 
classical limit e F r 3> 1, one obtains the solution Qq = 
of <S> of. However, as usual, this saddle-point is not 
unique. In particular, Q — TQqT^ 1 is also a solution, 
for any constant matrix T which is consistent with the 
fundamental symmetry (|l2|) for Q. Transverse fluctua- 
tions of Q away from the Q 2 — 1 manifold may be inte- 
grated out within the saddle-point approximation due to 
the large parameter vL d /r 3> 1, where L is the system 
size |2(J . Restricting attention to the manifold generated 
by the non- linear constraint Q 2 — 1, an effective action 
is obtained by allowing T to vary in space and expanding 
to second order in gradients of T, and first order in e, |A|, 
and 1/t s 



S[Q]=~ J drstv 



D(dQ) 2 -4i (e_o 
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(13) 



Here D = v^r/d represents the classical diffusion con- 
stant associated with the non-magnetic impurities. In 
particular, the quasi-particle DoS is obtained from the 
functional integral 



(i/(e,r)) ViS = -Re(str(of ®o p 



^ C Q(r))) ( 



(14) 



The numerical factor leads to a DoS of 4z/ for the system 
as |e| — > oo. This is because both the particle-hole struc- 
ture of the original Bogoliubov Hamiltonian and the spin 
each cause a doubling of the DoS. 

This completes the derivation of the intermediate en- 
ergy scale action. In the following section we will inves- 
tigate the action ( |T3| ) within mean-field theory and the 
influence of soft fluctuations on the low-energy properties 
of the system. 



D. AG Mean-Field Theory and Fluctuations 



The presence of symmetry-breaking terms in (|13|) orig- 
inating from the order parameter and the magnetic impu- 
rity potential means that a mean-field analysis is already 
non-trivial. To assimilate the effect of these terms, and 
to establish contact with the AG theory, it is necessary 
to explore the saddle-point equation. 



Varying the non-linear o-model action with respect to 
fluctuations of Q, subject to the non-linear constraint, 
one obtains the saddle-point or mean-field equation, 



Dd (QdQ) + i [Q, e_of <g> of + i\ A^™] 



6r s 



PH 
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0. 



With the Ansatz: 
Qmf = 0a 



of COsh 9 



1 sinh 9 



where the elements 9 = diag(#i, i9) BF are diagonal in 
the superspace, the saddle-point equation decouples into 
Boson-Boson and Fermion-Fermion sectors, and takes the 
form 



d 2 /( 9 + 2i [^cosh 9 - sinh 9 -( sinh(20) = , (15) 

a result reminiscent of the Usadel equation of quasi- 
classical superconductivity pj| ]. This is no coincidence: 
when subject to an inhomogeneous order parameter, the 
same effective action (O) describes the proximity ef- 
fect in a hybrid normal /superconducting compound [ p^[ . 
In the present context, the spatially homogeneous form 
of Eq. (15) should be combined with the self-consistent 
equation for the order parameter 

™9A \ - 
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Here <?a is the BCS coupling constant and 9 n indicates 
that the solution of the Usadel equation is taken at Mat- 
subara frequency e — > ie n . These two equations then 
coincide with the mean-field equations obtained by AG, 
and correctly reproduce the known form of E gap as spec- 
ified by Eq. (|). 

The AG solution is not unique: for e — > 0, the saddle- 
point equation admits an entire manifold of homoge- 
neous solutions parameterized by the transformations 
Q = TQupT' 1 where now T = 1 PH (g) 1 SP <g) t and 
t = 7(f~ 1 ) T 7 -1 . Soft fluctuations of the fields, which 
are controlled by a non-linear cr-model defined on the 
manifold T G OSp(2|2)/GL(l|l) (symmetry class D in 
the classification of Ref. pl|), control the low-energy, 
long-range properties of the gapless system giving rise 
to unusual localization and spectral properties. (For a 
comprehensive discussion of the physics of the gapless 
phase, we refer to Refs. [p2|-|2(| . ) The phenomenology of 
Class D in the context of the gapless phase of the present 
system will be the subject of a forthcoming paper [fj7| . 
Here we focus on the gapped phase, where Class D fluc- 
tuations are just one of a host of massive modes that are 
unimportant in the description of sub-gap states. 

This completes the formal description of the bulk su- 
perconducting phase. The solution of the AG mean-field 
equation provides an adequate description of the bulk 
extended states. Soft fluctuations around the AG mean- 
field describe phase coherence effects due to quantum 
interference. However, within the present scheme it is 
not yet clear how to accommodate sub-gap states in the 
gapped phase of the AG theory. To identify such states, 
it is necessary to return to the saddle-point equation ( |l5| ) 
and seek spatially inhomogeneous solutions. 



III. INSTANTONS AND SUB-GAP STATES 

Although the reduction and eventual destruction of 
the quasi-particle energy gap predicted by the AG mean- 
field theory can be reasonably justified on purely physical 
grounds, the integrity of the gap of the range < ( < 1 
is less credible. Once time-reversal is symmetry is broken 
and the protection of Anderson's theorem is lost, there 
remains no reason why a sharp gap should persist. Add 
to this the observation that the spin scattering rate must 
be subject to spatial fluctuations from the average value 
1 /t s , and one concludes that corrections to the DoS pre- 
dicted by the AG theory must lead to the appearance of 
sub-gap states analogous to "band tails" in a disordered 
semiconductor j[(J , [yJ . 

This analogy is of course not new [0j8| nor, as far as 
practical calculation in the present formulation is con- 
cerned, is it particularly deep. This is because all aver- 
ages have already been taken, so we can not look for an 
optimal fluctuation of some potential, as in the classic 
approaches to the study of band tail states in disordered 
semi-conductors 11 1. However, these studies hint at how 



one can proceed. 

Band tail states in semi-conductors can be studied 
within the same functional integral formulation. In par- 
ticular, the generating function of the single-particle 
Green function of a normal disordered conductor can be 
presented in the form of a supersymmctric field integral 



Z[0]= / £>(*,*) exp 



i I dr^ e + - — V(r) | * 

2m 



where, once again, the random impurity distribution is 
drawn from a Gaussian (^-correlated white-noise impurity 
potential. The optimal fluctuation method involves mini- 
mizing the action with respect to fluctuations in the fields 
"J and potential V. This involves seeking inhomogeneous 
solutions of the non-linear Schrodinger equation 



V(r) * = 0, 



2m 



where the corresponding optimal potential is determined 
self-consistently by the relation V(r) = — \^(y)\ 2 /2ttvt. 
In the supersymmetric formulation, band tail states are 
identified with supersymmetry broken inhomogeneous so- 
lutions of the saddle-point equation (see Cardy p8| and 
Affleck ^9|). Indeed, the anticipated exponential sup- 
pression of the DoS necessitates a breaking of supersym- 
metry to support a finite action. Here the phrase "su- 
persymmetry breaking" is potentially misleading. We 
use it only to refer to field configurations, ubiquitous in 
the problems under discussion here, that do not respect 
the parity between Bose and Fermi degrees of freedom. 
However, any such configuration is just one member of a 
degenerate manifold differing by supersymmetric trans- 
formations. The latter maintain the invariance of the 
generating functional Z[0] under global supersymmetric 
transformations. 

What does this tell us about the identification of opti- 
mal fluctuations and sub-gap states in the superconduc- 
tor? Following the analysis above, one might guess that 
sub-gap states are associated with inhomogeneous con- 
figurations of the $ field action. However, we anticipate 
that optimal solutions corresponding to sub-gap states 
are localized on a length scale in excess of the supercon- 
ducting coherence length. In the dirty limit, £ 3> £ ^> Xp, 
this implies that the localized sub-gap states are quasi- 
classical in nature. Their existence on the level of the \E' 
field action will be buried in the fast Xf oscillations of 
the wavefunction. To reveal the sub-gap states, we must 
first remove the fast short length scale fluctuations of the 
quasi-classical Green function and look for an equation 
of motion for the slowly varying envelope of the wave- 
function. But this is just the program of the usual quasi- 
classical method. 

The term "sub-gap states" is a little misleading in this 
context. Band tails are bound states of some rare poten- 
tial that sit by themselves below the bulk of the spec- 
trum. Each rare configuration that make the gap soft 
in the present case will give rise to many states beneath 
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the AG gap. Thus the term "gap fluctuation" , used in 
Rcf. jl5| to describe the zero dimensional SN system, may 
be more appropriate. 

As well as being quasi-classical in nature, the existence 
of sub-gap states is not affected by working in the dirty 
limit. As such, their existence must be accommodated 



in the non-linear cr-modcl functional (13) since the valid- 
ity of this description relied only on the quasi-classical 
parameter tpT ^> 1 and the dirty limit assumption. 
To identify sub-gap states in the present formalism, we 
should therefore investigate inhomogeneous solutions of 
the low-energy saddle-point equation in Q — the Us- 
adel equation p9[ . Such a solution should be thought 
of as defining an envelope for the quasi-classical sub-gap 
states. 

Therefore, let us revisit the mean-field equation and 
look for inhomogeneous solutions at energies e < E gap . 
To focus our discussion, let us begin by restricting at- 
tention to the quasi one-dimensional geometry. To stay 
firmly within the diffusive regime, we therefore impose 
the requirement that the system size L be much smaller 
than the localization length of the normal system £i oc ~ 
vLyjD, where L w denotes the cross-section. Later, in 



section [II C , we will generalize our discussion to encom- 



pass systems of higher dimension. Furthermore, since, 
over the interval < £ < 1, the quasi-particle energy 
gap varies more rapidly than the superconducting order 
parameter, we will neglect self-consistency of the order 
parameter. Taking self-consistency into account will not 
alter our qualitative findings, and will only weakly affect 
the quantitative results. 



A. Instantons in the Quasi One-dimensional 
Geometry 

To investigate inhomogeneous solutions of the mean- 
field equation (^) it is convenient to recast the equation 
in terms of its first integral 



where 



(d x/i 6f + V(9) = const, 



V{8) = U (sinh0 - — cosh 9 



(16) 



C cosh 29 



denotes the complex potential. Let us denote by 9 AG the 
values of 9\ and i9 at the conventional saddle point, and 
focus on an energy e below the gap predicted by the AG 
theory. Here Im 9 AG = ir/2 such that the mean- field 
DoS v AG (e) = 4^Rc cosh0 AG vanishes. The correspond- 
ing value of Re 9 AG depends sensitively on the energy, 
with Re 9 AG = for e = 0. 

Considering the Boson-Boson sector only, if we require 
that 9\{x — > ±oo) = #AGj what kind of inhomogeneous 
solution is possible? The values of 9\ at which d x 9i =0 
can be identified by considering the complex (dimension- 
less) potential function V(6\) from which we can deter- 
mine the endpoints of the 'motion' in the complex plane, 



just as one would use a real potential normally. By in- 
spection one may see that, on the line Im 9\ = tt/2, 
the potential is purely real. This is not the only con- 
tour where Im V = 0, but, by considering forces, it is 
not hard to see that either Im 9% — tt/2 always during 
the motion, or 9\ follows a trajectory with an endpoint 
at Im 9\ < 0. For reasons outlined below, we will dis- 
count this latter possibility. The former case amounts 
to considering "bounce" trajectories in the real potential 
V(m/2 + 0) = V R ((f>) where 



v R (<p) 



-4 ^cosh( 



— sinh0 +Ccosh20. (17) 



A typical potential is shown in Fig. 



v R - 5 



-10 



-15 




-2 2 4 



FIG. 2. Potential V K (<f>) = V(m/2 + 4>) for e/\A\ = 0.1 and 
£ = 0.2. The AG saddle point corresponds to the central max- 
imum. The saddle point marked <f>Q D is used in the analysis of 
the zero-dimensional problem (section IV). 



Now integration over the angles 9 is constrained to 
certain contours |^(J. Is the bounce solution accessible 
to both? As usual, the contour of integration over the 
Boson-Boson field 9\ includes the entire real axis, while 
for the Fermion-Fermion field, i9 runs along the imagi- 
nary axis from to Mr. With a smooth deformation of 
the integration contours, the AG saddle-point is accessi- 
ble to both the angles 9 p3| . By contrast, the bounce 
solution and the AG solution can be reached simultane- 
ously by a smooth deformation of the integration con- 
tour only for the Boson-Boson field 9\ (see Fig. ||). The 
bounce solution is therefore associated with a breaking of 
super -symmetry at the level of the saddle point. 
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FIG. 3. Integration contours for Boson-Boson and 
Fermion-Fermion fields in the complex 9 plane. The bounce 
solution for e — (labelled as 'b') is shown schematically. 

Thus we have identified an inhomogeneous saddle- 
point configuration for which the supersymmetry is bro- 
ken: 9\ executes a bounce whilst id remains at the mean- 
field value 9 AG . The symmetry broken solution then in- 
curs the (finite) real action 

S = 4 7 r^ w ( J D|A|) 1 / 2 5 (e/|AU) 
where, defining 4>' as the endpoint of the motion, 



Sft, = 



ci<h/K(0A G )-KW>). 



(18) 



0AG 



Now, as mentioned above, there exists a second pos- 
sibility for a bounce solution in which one moves away 
from 9 AG parallel to the imaginary axes. Indeed, such a 
solution would seem to be a natural candidate for the 
Fermion-Fermion field i9. However, since the endpoint 
for this trajectory lies at Re 9 < outside the integra- 
tion domain which runs from to n, this would seem to 
be excluded. 

As e approaches E gap from below, the potential (|l^ ) 
becomes more shallow, with the maximum merging with 
one of the minima when we reach the gap. Near the 
edge, up to an irrelevant constant, an expansion of the 
potential in powers of ((f) — 4> AG ) leads to the cubic form 



ap 



1/2 



+p(4> - </> AG ) 3 

where the dimensionless coefficients are specified by 

C-E'gap 



1/6 







1/3 



(19) 



(20) 



Note that, making use of Eq. (g), both of these coeffi- 
cients depend solely on the dimensionless parameter £. 
From this expansion, one can obtain an analytic solution 
for Scf,. To leading order in (-E gap — e)/|A| one finds 



S r h 



4 a 5 / 2 / E. 



15 (3 2 



\A\ 



5/4 



(21) 



This approximation is shown in Fig. ^ along with the ex- 
act result obtained by numerical integration. Note that 
the action vanishes exactly at the gap. For completeness 
we give the explicit form of the bounce solution 

1 

4>(x) - 4> A 



(3 cosh 2 (a;/2ro) ' 
where the extent of the instanton is set by 



r (e) 



*V2 V E, 



1/4 



(22) 



Indeed the size of the instanton is easily understood from 
the quadratic "stiffness" term in Eq. ( |l9| ) . Thus one finds 
that, while the overall scale is set by the superconducting 
coherence length £, the size of the droplet diverges both 
as e approaches E gap and, noting that a ~ (1 — C 2 / 3 ) 1 / 4 , 
as one approaches the gapless phase C — ► 1. 



1.5 





> 




\ A 

S 1 


-e- 


S3 1.0 

_o 


\ \ 


Act 


\ \ 
\ \ 




FIG. 4. Action S$ for C, = 0.5 obtained numerically (solid 
curve) together with the expansion in (-E gap — e)/|A| (dotted 
curve) as determined by Eq. (^9|). Note that the action van- 



ishes as e 
inset. 



The AG solution for the DoS is shown 



This completes the analysis of the saddle-point solu- 
tion together with the corresponding action. However, 
as this level we are presented with two problems: 

• the contribution of a second saddle point 
would seem to spoil the normalization condition 
(Z[0])y.s = 1, which should be preserved within 
the saddle point approximation; 

• confined to the line Im 9 = 7r/2, when substituted 
into the DoS source (fl4|), the bounce configuration 
does not appear to generate states! 

The resolution of both problems lies in the nature of the 
fluctuations around the symmetry broken mean-field so- 
lution. These field fluctuations can be separated into "ra- 
dial" and "angular" contributions. The former involve 
fluctuations of the diagonal elements 9, while the latter 
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describe rotations including those Grassmann transfor- 
mations which mix the bf sector |3^]. Both classes of 
fluctuations play a crucial role. 



Defining 



C A 

b rc^T- 1 



B. Fluctuations 

Before turning to the technical analysis, let us out- 
line qualitatively the influence of the fluctuations around 
the mean-field. As usual, associated with radial fluc- 
tuations around the bounce, there exists a zero mode 
(due to translational invariance of the solution), and a 
negative energy mode. The latter, which necessitates 
a 7r/2 rotation of the corresponding integration contour 
to follow the line of steepest descent (c.f. Ref. (3|] and 
see below), has two effects: firstly it ensures that the 
non-perturbative contributions to the local DoS are non- 
vanishing, and secondly, that they are positive. Turning 
to the angular fluctuations, the breaking of supersymme- 
try is accompanied by the appearance of a Grassmann 
zero mode separated by a gap from higher excitations 
which restores the global supersymmetry (c.f. spin sym- 
metry breaking in a ferromagnct of finite extent). The 
zero mode ensures that the symmetry broken inhomoge- 
neous saddle-point configurations respect the normaliza- 
tion condition (Z[0])y.s = 1. 

To formally investigate the fluctuation determinant, let 
us implement the rational parameterization 



1 + iP 
1 -iP 



R~ 



(23) 



where the condition P = o*^YP T T 1 a™ with 



cc 

03 7 



0i°° 



is imposed b y (P-2|) . With this choice, a variant of that 
used in Ref. |32] , the measure is trivial [EOj . R rotates Q 



from the metallic saddle point a™ 
configuration: 

R(x) = exp 



crS c to the bounce 



6{x) 



where A = -TA T T 
Q 2 = 1 requires [a 



and B = —TB T T l , the condition 
P]+ = 0, or 



'3 ! 



Thus the matrix field C is off-diagonal in the CC space, 
while A and B are diagonal. In fact, the field fluctu- 
ations C describe the low-energy quantum interference 
effects — the soft Class D modes. However, the fluctu- 
ations contained within the fields C are oblivious to the 
supersymmetry breaking. We will therefore deal only 
with the fluctuations that are parameterized by A and 
B. (Moreover, we will neglect the massive spin triplet 
fluctuations.) 

With the explicit parameterization [B3| 



.4 



B 



Sb 





')A 













VA 


-Va 





Sp 








VA 





-Si 


— s* 





Vb 













Vb 


-Vb 













VB 
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an expansion of the action ( |l3| ) to quadratic order around 
the bounce configuration 



6{x) 



h(x) 




obtains 



5 = 4^i^(i^|A|) 1 / 2 (^ + ^ q ) 



where represents the contribution from the saddle- 
point alone (lq), and 



du 



{\ds B \ 2 + |3 Sf | 2 ) + F'(0)4 2 + v"{4>K 2 + V'( 



, 2 



V" 



0« 



r/2 



(24) 



where s b/f ee s' b/f + is% /p , C± = {vb ± Va)/V2, £± = 

(fjB ± 77a)/v2, u = x/£ and 9 AG = z7r/2 + <p AG . Here the 
various potentials are given by 



V'(4>) = 2 cosh(/> 



e 

IAI 



sinht/) — Ccosh2</> 



r 



V"{4>) = (d(f)) 2 + 2 ( cosh <p - -— sinh <f> - £ sinh 2 4> 



and 



(dej)) 2 

V± (^) = — H cosh + cosh <j) AG 







- — (sinh 



3 ) — — (cosh 20 + cosh20 A 



sinh y-*ag 

C (sinh sinh 0ag T cosh cosh </>ag) 



As a check, let us consider the conventional saddle 
point. At e = (for simplicity), (f> = cf> AG = 0, and 
the quadratic action assumes the form 



S q = du (\ds B \' + \8s F \ z + d£ + dt+ + d£-d£- 



+2 



du 



s B 

l 2 



/2 



This action is manifestly supersymmetric and, therefore, 
performing the integrations over all fields gives unity. 
(Moreover, since £ < 1 in the gapped phase, the inte- 
gral is manifestly convergent.) This in why the usual AG 
DoS is just given by evaluating the source ([l4]) at the 
supersymmetric AG saddle point. 

Since V'((f> AG ) and V"((f> AG ) are both positive definite, 
integration over the Fermion-Fermion degrees of freedom 
merely generates some (weakly energy dependent) posi- 
tive prefactor. The Boson-Boson sector is more interest- 
ing. In particular, it is straightforward to verify that the 
action for the components s' B simply reflects the longitu- 
dinal variation 



\ 1 < U 



5 2 S 



with 2s' B — A6*i. Thus the action for s' B is the one that 
could have been written down from the outset: it is that 
of the fluctuations of 9\ around the bounce instanton dis- 
cussed in section III A. To address the influence of this 



class of fluctuations we can draw on the standard litera- 
ture [glj. 

To perform the Gaussian integration over the fields s' B , 
we form the expansion 

s 'b( x ) = ^2a n y> n (x) 

n 

in terms of the eigenfunctions ip n of the quadratic opera- 
tor in the action for s' B . Now the action for s' B exhibits a 
zero mode ifi ~ d6\ due to translational invariance of the 
action (the bounce can be positioned anywhere in space) . 
This zero mode must be accommodated by the introduc- 
tion of a collective coordinate which in turn introduces 
a Jacobian factor associated with the change of variables 
from a\ to the collective coordinate. Furthermore, since 
the instanton is a bounce, the zero mode (f \ has a node, 
and hence there exists a ground state of V with nega- 
tive energy. This requires the contour for the integration 
variable clq to be deformed away from the real axis - 
it takes a right turn at zero and heads in the negative 
imaginary direction. In the context of the contour drawn 
in Fig. ||, this is because we have to get back to the real 
axis for 6\. 



This deformation of the integration contour has a pro- 
found consequence. While the contribution to the sub- 
gap states from the bounce solution alone vanishes (recall 
that the instanton was confined to the line Im 9\ = tt/2, 
so that cosh 0i remains imaginary), the rotation of the 
integration contour introduces a factor of i resulting in 
an imaginary contribution to the Green function below 
the gap. This is turn would signify a non-zero sub-gap 
DoS. The mechanism operates in the related context of 
Landau levels broadened by disorder and discussed by 
Efetov and Marikhin |34| (see also an earlier paper by 
Affleck H). 

It is reassuring to note that the deformation of the in- 
tegration contour, which is constrained to return to the 
undeformed Bosonic contour, is unambiguous, and gives 
rise to only positive definite contributions to the DoS. 
Finally, as usual, since the contour only runs over half a 
Gaussian for ag, its contribution yields a factor of 1/2. 




FIG. 5. Spatial dependence of the Grassmann zero mode 
(Pq for e = and £ = 0.5 together with the potential V- 
(scaled by |A|) that binds it. 

This completes the discussion of the contribution to the 
functional integral from the Bosonic degrees of freedom. 
Finally we turn to the role of the Grassmann fluctua- 
tions. Considering the supersymmetric structure of the 
action, it is immediately clear that supersymmetry break- 
ing must be accompanied by the existence of zero mode 
in the Grassmann sector which restores the supersymme- 
try. Indeed, an inspection of the potential V- identifies 
a zero energy eigenfunction gapped from the others (see 
Fig. 1)11]. 

This situation may be compared with the analysis of 
Andreev and Altshuler |pq| , who identified a stationary 
phase saddle-point that determines the oscillatory part 
of energy level correlation functions for normal diffusive 
conductors. Again this saddle-point is more correctly a 
manifold of points related by supersymmetry, but the fact 
that it is spatially uniform guarantees that the spectrum 
of Grassmann fluctuations is truly gapless. 
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This zero mode indicates the existence of a degener- 
ate manifold parameterized by a Grassmann coordinate. 
The integration of unity over this Grassmann coordi- 
nate yields zero, which ensures the necessary normaliza- 
tion condition (Z[0])v,s = 1 is met. A non-zero result 
can be obtained only if there is a prefactor (or source 
term) that breaks supersymmctry. Making use of the 
DoS source ( |l4]) , we expand in P to find the lowest order 
part quadratic in the Grassmanns (at order P 2 ) 

I 



v.s 



4z/ 



str (af ® cr™ <g> afQ{x)) ~ 

%i{svnh(j>(x) - sinh0 AG ) x (£ + (x)£+(x) + |_ (x)£- (x)) . 

From this result, we see that the local sub-gap DoS re- 
mains non-zero only in the vicinity of the instanton where 
the supersymmetry is broken. 

Thus, taking into account Gaussian fluctuations and 
zero modes, one obtains the non-perturbative, one in- 
stanton contribution to the sub-gap DoS: 



(— i\K\) / dx i(smh(f>(x) — sinh^> AG ) | y? (a;) | ' 



LSrh 



exp 



-4irvL w ^D\A\S 4 



(25) 



where the factor y/LS^Jt; represents the Jacobian as- 
sociated with the introduction of the collective coordi- 
nate |n], — i\K\ is the overall factor arising from the 
non-zero modes, and the Grassmann zero mode wave- 
function ifg is normalized such that 



dx\ip Q | 2 = 1. 



Here we have assumed that the integration only con- 
tributes to the positive prefactor. We have checked this 
numerically for a few cases and expect a general state- 
ment could be made by moving to a different parameter- 
ization. 

Eq. ( p5| ) is the main result of this section. Note the 
non-perturbative nature of the result, both in the cou- 
pling constant g~ x of the cr-model, and the (dimension- 
less) spin scattering rate £. 



I 

that the DoS depends exponentially on the parameter 
<±Kg{ilL) d - 2 S^ where 



'TP _ c \(S- d )/ i 
S, = a d C 2/3 (l C 2/ T (2+d)/8 ( ^ft 1 ) • (26) 



Here g = vDL d 2 denotes the bare dimensionless con- 
ductance of the normal system, and a<j is a numerical 



constant (ai = 8 V24/5) In particular, the exponent 
depends linearly on the energy separation from the gap 
in two dimensions. 



C. Sub-gap States in Dimensions d > 1 

The calculation above was tailored to the considera- 
tion of the quasi one-dimensional geometry. The gener- 
alization to higher dimensions follows straightforwardly. 
In particular, it is necessary to seek inhomogeneous solu- 
tions of the saddle-point equation (|l6| ) where the gradient 
operator must be interpreted as the higher dimensional 
generalization. Generally, this equation must be solved 
numerically. However, for energies e in the vicinity of the 
gap -Egap, an analytic expression for the energy scaling 
can be obtained. 

Using the approximation to V K [<f>] ( |l9| ) valid when 
(E-gap — e)/|A| -C 1, the exponential dependence of the 
sub-gap DoS can be deduced in higher dimension. In 
this limit, dimensional analysis of the cubic equation of 
motion yields the scaling form 

^r)-^ AG (e) = |/(r/r ), 



where ro is the characteristic length defined by Eq. (22) 



D. Numerics 



When substituted back into the action, one finds 



To assess the validity of the approximations used in 
obtaining the results above we have investigated numeri- 
cally the DoS in the vicinity of the mean-field energy gap 
using a numerical diagonalisation of a non-interacting 
tight-binding Bogoluibov Hamiltonian with disorder in 
both the on-site potential matrix elements and in the spin 
impurity scattering potential. Taking a two-dimensional 
lattice of size 22 x 22 site with an on-site disorder taken 
from the range of e e [— 3,3] measured in units of the 
hopping matrix element, |A| = 1 measured in the same 
units, and various strengths of the spin impurity poten- 
tial, the DoS is shown in Fig. ^. Notice that, as the 
strength of the spin impurity potential is increased, the 
quasi-particle energy gap is quenched. 
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FIG. 6. Variation of the DoS for a weakly disordered 
two-dimensional tight-binding Bogoluibov Hamiltonian in the 
presence of magnetic impurities. The data is shown for in- 
creasing strengths of the magnetic impurity potential. At the 
lowest energies the data points are shown as open circles to 
emphasize the fine structure. Notice that in the phase where 
there is predicted to be a gap in the mean-field theory the 
DoS shows a tail extending in the sub-gap region. When the 
quasi-particle energy gap is fully suppressed, the DoS shows 
an upturn at very low energies which is compatible with the 
renormalisation due to quantum interference processes pre- 
dicted by the class D theory. 

For the weakest magnetic impurity potential, we have 
expanded the region in the vicinity of the mean-field en- 
ergy gap. Fig. shows a exponential scaling of the 'sub- 
gap' DoS with an exponent which depends linearly on 
energy as predicted by the theory above. 




0.2 0.3 

e/[A| 



0.4 



FIG. 7. Variation of the DoS in the vicinity of the 
mean-field energy gap. The data is shown as open circles 
and is compared to an exponential fit of the data shown as a 
solid curve. Notice that the DoS depends exponentially on the 
energy difference as predicted by the two-dimensional version 
of the theory. 



IV. ZERO DIMENSIONAL PROBLEMS AND 
UNIVERSALITY 

In the previous section we considered the instanton 
contribution to the sub-gap DoS in the infinite system. 
For completeness let us now consider the zero dimen- 
sional case that obtains when ro, the size of the instan- 
ton, exceeds the system size L, which will happen when e 
approaches close enough to -E gap from below in any finite 
system. In this limit one can clearly not fit an instanton 
inside the system. Leaving aside the practical relevance 
of this situation, theoretical motivation is provided by a 
recent paper || that explored this regime using a random 
matrix analysis. 

However, before turning to the consideration of the 
zero-dimensional limit of the present problem let us first 
try to draw some intuition from a closely related in- 
vestigation of a different system comprised of a normal 
quantum dot contacted to a superconducting terminal 
(as shown in Fig. ||). In such a geometry it is well es- 
tablished (see e.g. Ref. J37j]) that the proximity effect 
induces a gap in the DoS of the normal dot. Indeed, in 
Ref. [[57J , the integrity of the gap is proposed as a signa- 
ture of irregular or chaotic dynamics inside the dot. Now 
near the gap edge the DoS of the dot takes the singular 
form 



v{e > E g 




(27) 



with E gap = cNS and A g = c'N^-^S, where c « 0.048 
and c' « 0.068, 8 = 1/L d u denotes the single particle 
level spacing, and N is the number of fully transmitting 
channels in the contact. 

However, the location of the gap edge relies on a mean- 
field analysis of the coupled system. In Ref. Jl5| Vavilov 
et al. have argued that optimal fluctuations of the impu- 
rity potential give rise to gap fluctuations. The hypoth- 
esis introduced in Ref. |l5| is that the spectral statistics 
near a gap edge are universal. This allows a random ma- 
trix theory analysis of gap fluctuations and leads to the 
following expression for the sub-gap DoS, 



v{e < E gap ) 



exp 



3/2' 



(28) 



Now the AG mean-field solution for a superconducting 
quantum dot with magnetic impurities also predicts the 
existence of a square root edge (see Eq. (|J) and inset of 
Fig. ||). Then, when recast in the form of Eq. (|27|), it 
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is pertinent to ask whether the expression for the sub- 
gap DoS coincides with Eq. ( ^8|) in the zero dimensional 
limit. This is the situation addressed in Ref. j8j. 

In the following we will show that the universal re- 
sult (^||) is explicitly recovered by the present theory. 
Moreover, in doing so, we will expose the origin of the 
universal structure reported in Ref. |l5| and describe its 
implications for universality of the d > problem. 



part of the contour that gives a finite sub-gap DoS. In 
the zero-dimensional case we are spared having to think 
about the problem in function space. The action is pro- 
portional to the potential of Fig. ||. The correct contour 
thus passes through the maximum of the potential (mini- 
mum of the action) corresponding to the usual AG saddle 
point, and turns away from the real <f> axes (i.e. the line 
Im 9\ — 7r/2) at the minimum of the potential 




FIG. 
leads 



Metallic quantum dot coupled to superconducting 



A. Superconducting dot with magnetic impurities 

Let us therefore consider explicitly the action of a su- 
perconducting grain in the presence of a weak magnetic 
impurity potential. When tq 3> L only the zero spatial 
mode contributes significantly to the action ([l3]). In this 
limit, the action assumes the zero dimensional form 



5[Q] = -str[( e _a 3 cc 



°3 Ql 



+ 



" str (Qct™ ® ct sp ) 2 



24tJ 



(29) 



where, as usual, S denotes the single-particle level- 
spacing. As in the higher dimensional problem, a varia- 
tion of the action with respect to Q obtains a mean-field 
equation, now without spatial variati on. P arameterizing 
the saddle-point equation as in section II D , we obtain the 
zero-dimensional Usadel or AG mean-field equation fllf 



2i ^cosh 9 



sinh0 - Csinh(26») = 0. 



From this equation, we can identify the usual AG solu- 
tion which in turn recovers the AG phenomenology. 

The inclusion of bounce configurations in the previ- 
ous calculations was based upon the observation that, 
although the contribution they make is exponentially 
small, they are the least action configurations on the 




0Cto(e) = <?W(e) 



(marked in Fig. ||). This part of the contour, parallel to 
the imaginary axes, gives a contribution to the DoS, and 
the second saddle point is in fact a maximum on this por- 
tion by analy ticity. Following the same arguments as in 
section III A , this solution is inaccessible to the Fermionic 
contour. Wc find that the sub-gap DoS near -E gap in the 
zero-dimensional case scales as 



v(e < £ g ap) 



exp 



4tt |A| a 3 



ap 



3/2' 



where a and j3 are the coefficients defined in Eq. (|2C 
We note that the general result for the energ y dep en- 
dence of the exponent written down in section III C for 
dimensions d > 1 applies also for d = 0. 

To establish contact with the universal result given in 
Eq. (^) it is helpful to recast the result in a modified 
form. To do this we note that, in the vicinity of the 
mean-field gap edge, the DoS can be expanded as III] 



v AG (e > ff ga p) 
42/ 



£- 2 /3(l - £2/3)-l/4 



E, 



1/2 



Then, if we define 




£2/3n-1/4 



the mean-field DoS can be brought to the form of 
Eq. (p7|) , and the sub-gap DoS takes the universal 
form pal 



V(e < £gap) 



exp 



E, 



gap 



3/2' 



(30) 



Leaving aside a spurious (yet systematic — see below) 
factor of 2 p9[ |, the sub-gap DoS obtained above coin- 
cides with the universal expression shown in Eq. (pq). 

Once again, to assess the validity of the approxima- 
tion scheme that leads to the universal result above, we 
have looked for the scaling in a numerical investigation of 
a random matrix Hamiltonian with the same symmetry. 
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In the sub-gap region, Fig. ||| shows a good fit of the DoS 
to the predicted exponential scaling. 

The rescaling of the DoS above and the appearance of 
the universal form suggests that we should revisit the d- 
dimensional result and look for a similar rescaling. From 
Eq. ( p6| ) it is straightforward to verify that in this case 



v(e < £ gap ) 



exp 



-ad 



ro(e) 



ap 



3/2' 



where represents some numerical constant, and r is 
the characteristic length defined by Eq. (^) . Finally, by 
defining 



-3/2 



5(e) 



_ A -s/a 



where 5(e) — l/(vr^(e)) is the level spacing inside a re- 
gion of size ro, the volume dependent prefactor can be 
absorbed into the expression and DoS can be brought to 
the form 



v(e < £ gap ) 



exp 



-ad 



3/2' 



revealing a simple relation between the d = and d > 
problems. 

The coincidence of Eqs. (|2^) and ([30|) indeed suggests 
that gap fluctuations are universal. To expose the origin 
of the universal scaling within the present formalism, let 
us consider the quantum dot geometry of Fig. ^| within 
the (7-model scheme. 
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FIG. 9. DoS for a random matrix Bogoluibov Hamiltonian 
with a magnetic impurity component. Data is shown for three 
different ensembles where only the rank of the matrix is varied 
(M = 50, 100 and 200). The DoS has been rescaled by the 
average level spacing so that, above the mean-field gap edge, 
the data collapse. As predicted, with this rescaling, the width 
of the tail region decreases with increasing size of the random 
matrix. The data set for M — 200 is fitted to an exponential 
with a 3/2 power. The corresponding fit is then made with 
no adjustable parameters to the other two data sets. 



B. Quantum dot contacted to a superconductor 



Following Ref. [pL 5| , let us consider a normal quantum 
dot contacted to a bulk superconducting terminal with 
order parameter A. Taking N fully open channels be- 
tween the dot and the lead, it is straightforward to show 
that, in the zero-dimensional limit, the effective action of 
the hybrid SN system takes the form 



S[Q] 



N 

—str [ln(l + QlQ)] 



(31) 



where Ql represents the Q-matrix of the superconduct- 
ing terminal (see below), and Q represents the zero- 
dimensional supermatrix field within the normal region. 
The term containing the logarithm describes the coupling 
to the lead — see Ref. |g(J f° r a derivation. Even for a 
diffusive dot, we note that retaining the full logarithm 
is the essential element in the correct treatment of the 
zero-dimensional limit that obtains when D/L 2 3> N6, 
a point to which we will return later. For a clean super- 
conducting lead, at energies e much less than the bulk 
superconducting gap A, one can set Ql = <?i H - 

As usual, to obtain the mean-field expression for the 
DoS it is necessary to minimize the action with respect to 
variations in Q. Doing so, one obtains the saddle-point 
equation 



f[Q,(i 



QlQ)- 1 Ql} = 



As usual, applying the Ansatz that the saddle-point so- 
lution is contained within the diagonal parameterization, 

Q = 03° O cr™ cosh 9 + ia\ n sinh 9, 

the saddle-point equation takes the form 



7re . - N cosh# 

— sinh 9 + — 

5 2 1 + i s i n h 9 



0. 



(32) 



A straightforward analysis of the symmetric saddle- 
point solution leads to the mean-field result for the 
DoS shown in Eq. ( p7| ) (remembering that now v(e) = 
2v Re cosh# MF (e) as there are no spin degrees of freedom 
in this problem). In particular, one can straightforwardly 
determine E gap by setting cosh 9 MF to be imaginary. Thus 



sinh 9x 



-ib for real b and (p2|) gives 



e(b) 



1 



N51 b 
~2^b\ 6 + 1 



The extremum of this function gives the largest en- 
ergy corresponding to a real value of b. This occurs at 
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b = (1 + v5)/2 = 1 + 7, where 7 is the golden mean, and 
yields E gap = (N6/2ir)j 5 / 2 » 0.048^(5 as required. 

As with the case of magnetic impurities, to explore 
the influence of gap fluctuations, it is necessary to seek 
the symmetry broken saddle-point configuration. How- 
ever, for reasons outlined below, it is possible to identify 
a pattern which implies the universality of the resulting 
analysis: 




It is this symmetry broken saddle-point which controls 
the sub-gap DoS and leads to the universal scaling form 
proposed in Ref. fLSfl , This generalizes the arguments 
applied to the superconducting dot with magnetic impu- 
rities. 



C. Universalities 

At first glance the situations considered in the two pro- 
ceeding sections are rather different. The actions ( p9| ) 
and (plj) would seem not to have much in common. How- 
ever, a simple and general argument may be established 
to reveal the universal character. As before, defining 
^mf(e) = i7r/2+^> M p(e)i the mean-field DoS for the SN de- 
vice is given by ^ MF (e) = 2v Im sinh0 MF (e), where <fi 
determined by the condition SS/5<fi[(f> = MF ] = 0. Since 
the DoS displays a square root singularity described by 
Eq. (p7|), the (saddle-point) action near the edge is con- 
strained to be of the form 



s[4>] 



-k str 



1 



E v 



A3 



where s(e) = sinh^(e) — sinh</> MF (£ , g ap ). Here, the ele- 
ments <p = diag(</> BB , 4>ff) and s = diag(s BB , s FF ) are diag- 
onal in the superspace. (As one may check, a variation of 
the action for e > -E gap obtains the symmetric mean-field 
solution 



E, 



gap 



2tt 



AS 



which in turn recovers the expression ( p7| ) for v(e).) 
Moreover, since the term containing s is linear in the 
energy, we can determine the value of k from the knowl- 
edge that e appears in the action as (2n£+/S) sinh</>. (It 
is this term that can more generally contain the Dyson 
index '/3', which therefore appears in the general expres- 
sion for gap fluctuations described in Ref. |15|.) In the 
present case, we thus have k = (27rA g /<5) 3 . 

Now, as discussed in the previous section, when e < 
-Egap there exists two saddle-point solutions at 




As before, one of these solutions (s_(e) ~» </> MF (e)) is asso- 
ciated with the conventional symmetric mean-field solu- 
tion while the other represents a second saddle-point ac- 
cessible only to the Bosonic contour. Taking this second, 
symmetry broken saddle-point into account (i.e. setting 
,s BB = s + and s FF = s_), one obtains the saddle-point 
action 



D. Discussion 



Following on from this discussion, to conclude this sec- 
tion, let us make two remarks which bare on the univer- 
sality of the general scheme. The first of these remarks 
concerns the integrity of the scaling of the sub-gap DoS 
when different impurity distributions are taken into ac- 
count. The second remark concerns the extension of the 
ideas above to the consideration of the hybrid SN system 
beyond the zero-dimensional regime. 

Firstly, for the superconductor with magnetic impu- 
rities, one can generalize the arguments above to show 
that the energy scaling of the sub-gap DoS even in the 
d-dimensional case is insensitive to the nature of the ran- 
dom impurity distribution. This is in contrast to Lifshitz 
band tail states in semi-conductors where the energy scal- 
ing depends sensitively on this distribution. To under- 
stand this, let us suppose that the distribution of mag- 
netic impurities JS(r) is not Gaussian (5-correlated, as we 
assumed throughout, but obeys some arbitrary statistics 
defined by a probability functional P[JS(r)]. When the 
ensemble average over JS(r) is performed one would ob- 
tain in the ^-field action a contribution of the form 



In 



exp 



dr *JS ■ cr sp * 



C[*a sp *(r)] , 



which defines C[---], the generating functional of con- 
nected correlators of JS(r). Though this is in general 
a very complicated and indeed non-local functional of 
^aty(r), one can in principle find a local Q-ficld action by 
including pairings only at coincident points, justified by 
the assumption (i/£,) d <C 1 about the non-magnetic dis- 
order. The mean-field description of this system then fol- 
lows from the homogeneous solution of the saddle-point 
equation, an Usadel equation like ( |l5|) with some poten- 
tial. Generally this potential will have the same char- 
acteristics as the real potential of ([lj]) plotted in Fig. || 
on the line Im 9\ = n/2. The central maximum is due 
to the |A| term; the upturn at large <j> arises from the 
small pair-breaking part, and the asymmetry comes from 
the e term. Now, if mean-field theory leads to a square- 
root singularity in the DoS (a circumstance which can be 
avoided only by a special tuning of parameters), one can 
expect that increasing the energy leads to the maximum 
merging with one of the minima according to 
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Vr[4>] 



1/2 



(j) MF ) 2 + P(<p - (f> MF ) 3 



with a an d 8 chosen appropriately. Then the analysis 
of section III applies. In particular the scaling of the 
exponent with ((E gap — e)/|A|)( 6 ~ d )/ 4 is expected to be 
universal and independent of the details of the magnetic 
impurity potential. 

Now let us turn to the generality of the present scheme 
in describing 'gap fluctuations' in extended hybrid su- 
perconductor/normal systems. The latter has been dis- 
cussed in a very recent paper by Ostrovsky et al. flfj| . In 
this work, the authors developed an instanton approach 
analogous to that employed in the magnetic impurity 
system here to estimate the profile of gap fluctuations 
in the d-dimensional SNS system. Now from the dis- 
cussion above, it is possible to expose the relation be- 
tween these two works: in the SNS system, the energy 
gap induced in the normal region due to the proximity 
effect is determined by the Thouless energy defined as 
Et ~ 1/Tdweii, where Tdwcii is the time required for elec- 
trons in the normal region to feel the presence of the 
superconductor J37j . The Thouless energy is determined 
by Et ~ mm{D/L 2 , TN5} where T is the transparency 
of the contact to the superconductor (T = 1 in the anal- 
ysis of the zero-dimensional system above). 

In the diffusive limit D/L 2 <C TN5, at the mean- 
field level, the position of the quasi-particle energy gap is 
found by solving the Usadel equation with the appropri- 
ate boundary conditions JO]. As a result one obtains a 
square root singularity in the DoS. In this case the mean- 
field solution is itself inhomogeneous. The sub-gap cor- 
rection is found by identifying a second inhomogeneous 
instanton configuration that breaks supersymmetry at 
the level of the action fjoj . Both solutions merge at the 
mean-field gap. Now, following the arguments above, it 
is simple to see how the phenomenology of Ref. |4(J fits 
into the same general scheme: in this case the relevant 
coordinate of Q interpolates between the inhomogeneous 
mean-field solution and the instanton. The result is a 



sub-gap DoS which assumes the familiar form of Eq. (30), 
with appropriately defined geometry dependent param- 
eters -Egap and A g . Naturally the introduction of dx 
transverse dimensions gives the expected energy depen- 
dence of (£gap — eY 6 ~ d± ^ 4 in the exponent. 

In the opposite limit D/L 2 3> TN5 (not considered 
in Ref. pO{), gradients of Q are heavily penalized and 
the coupling to the lead must be retained in its 'log- 
arithmic form', with Q being taken as constant in the 
dot. (Indeed, the logarithm is crucial to reproduce even 
the mean-field expression for the DoS ( p7| ) with the cor- 
rect coefficients.) This is the true zero-dimensional limit 
treated above. As we have seen, with this action, one re- 
covers the known universal expression for the spectrum 
of gap fluctuations below the mean-field edge. 



V. CONCLUSIONS 



To conclude, we have developed a quasi-classical field 
theory describing a superconductor in the dirty limit 
with weak magnetic impurity scattering. The Abrikosov- 
Gor'kov mean-field treatment of this system, showing a 
diminished but hard gap in the DoS, can be straightfor- 
wardly recovered as a homogeneous saddle-point of the 
effective action. Where zero DoS is predicted by the 
mean-field theory, there exist spatially inhomogeneous 
saddle-point configurations that break supersymmetry at 
the level of the action. A careful analysis of fluctuations 
around these instanton configurations demonstrate how 
supersymmetry is restored by a manifold of equivalent 
configurations parametrized by a Grassmann coordinate, 
but more importantly how the configurations give rise 
to a finite, though exponentially small, DoS. In contrast 
to band tail states in semi-conductors, the quasi-particle 
nature of the sub-gap states leads to universality of their 
properties. 

Finally, let us remark on the connection of the re- 
sults presented above to related problems in the liter- 
ature. The resulting expression for the DoS (^6|) was 
found to be non-perturbative in the a-model coupling 
l/<7, which measures the strength of non-magnetic disor- 
der. We note that other non-perturbative results in dis- 
ordered systems have been obtained by related instanton 
calculations. As well as the investigation of tail states in 
semi-conductors p8| , a supersymmetric field theory was 
developed by Affleck [^9| (see also Refs. |Q) to inves- 
tigate tail states in the lowest Landau level. There it 
was shown that tail states correspond to instanton con- 
figurations of the ^ -field action (c.f. Ref. |2£|). It is also 
interesting to compare the present scheme with the study 
of 'anomalously localized states' Q (see also, Ref. Q)- 
There one finds that long-time current relaxation in a 
disordered wire is also associated with instanton config- 
urations of the cr-model action. Finally, a Lifshitz argu- 
ment has been applied on the level of the Usadel equation 
in the study of gap fluctuations due to inhomogeneities 
of the BCS interaction Q. 

Although we have focussed largely on the question of 
tail states below the mean-field quasi-particle gap of a 
superconductor with magnetic impurities, we expect the 
instanton approach developed here to be more widely ap- 
plicable. Indeed, in this work we have shown the inti- 
mate connection between the study of sub-gap states in 
the magnetic impurity problem and gap fluctuations in 
hybrid superconductor/normal structures. Furthermore, 
the same instanton approach describes gap fluctuations 
in superconductors with a quenched inhomogeneous dis- 
tribution of the BCS coupling constant ^3|, as well as 
quasi-two dimensional superconducting films subject to 
strong in- plane magnetic fields [Q . In both cases the lat- 
ter are described by a mean-field theory which assumes 
the Abrikosov-Gor'kov form. 

More speculatively, it seems likely that the same gen- 
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eral scheme can be employed to study the influence of 
optimal fluctuations on the nature of bulk transitions. 
For example, in the present system, the transition to 
bulk superconductivity with reducing magnetic impu- 
rity concentration will be preempted by the nucleation 
of superconducting islands or droplets within the metal- 
lic/insulating phase (c.f. Ref. (4{|). Similarly, the Stoner 
transition to a bulk itinerant ferromagnet in a disordered 
system will be mediated by the formation of a droplet 
phase in which islands become ferromagnetic |46f| . In 
both cases, we expect these 'droplet phases' to be associ- 
ated with inhomogeneous instanton configurations of the 
corresponding low-energy action. 
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